Quantum communication through Jaynes-Cummings-Hubbard 

arrays 

R. Chakrabarti^ and G. Sreekumari^^^ 

"'^Department of Theoretical Physics, University of Madras, Guindy Campus, Chennai 600 025, India 
^Department of Physics, Loyola College, Chennai 600 034, India. 



^ PACS numbers: 42.50Pq, 42.50Ex, 03.67.Hk, 02.30. Gp 

CN ', Abstract 

?— j ' We study the dynamics of an one dimensional array of Jaynes-Cummings-Hubbard system of arbitrary 

number of coupled cavities, each containing a two level atom that interacts with a field mode. In 
particular, we consider propagation of a single excitation quantum state for two different couplings of 
the photonic modes of the adjacent cavities, namely, a translation invariant closed chain of uniformly 
coupled cavities, and also a linear chain with nonuniform parabolically varying intercavity coupling where 
the interaction Hamiltonian is associated with the Jacobi matrix of the Krawtchouk polynomials. Using 
^ a description via the delocalized atomic and field modes we observe that for a large detuning of these 

CIh' two degrees of freedom atomic excitations propagate without populating the field modes, and vice versa. 

For the near-resonance scenario between these modes the atomic excitations, say, while propagating mix 
P3 , with the photonic states. In the context of the parabolic coupling between photons of adjacent cavities 

an arbitrary element of the time-dependent correlation function between two arbitrary cavities may be 
^ I expressed in closed form for dominant values of the detuning parameter, when an exact transmission of 

the quantum state at pre-specified times is realized. 
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> ! I Introduction 

I Recently considerable theoretical and experimental attention has been devoted to a class of models involving 
coupled optical cavities where photons interact with embedded two level atoms. On the experimental side 
these models have been facilitated by experimental advance in photonic crystals [I] , optical microcavities [5] , 
■ and superconducting devices [3]. These coupled cavity structures have potential applications as quantum 
optical simulators [H [5] of condensed-matter phenomena making it feasible to study many-body effects 
such as phase transitions where the particles of interest are photons. The coupling between cavities also 
provides a setting for the system being considered as a wave guide in the context of distributed quantum 
information processing [S]. Allowing for control and measurement at individual lattice sites proposals have 
^ I been put forward for generating entangled photonic states [7] , creation of cluster states [51 [3] , and transfer of 
H ■ a quantum state along an array of polaritonic qubits [TO] . Towards investigating the transmission of quantum 
- - -' states in an one-dimensional coupled array of cavities much study has been done in the context of Jaynes- 
Cummings-Hubbard model [II) . The atom-cavity photon interaction is described by the well- understood 
Jaynes-Cummings model [T^] that relies on rotating wave approximation. Tunneling of photons between 
adjacent cavities is introduced via a hopping parameter. As the number of excitations in the system remain 
conserved it is, from the point of view of propagation of quantum states, of importance to develop a detailed 
understanding of the time evolution of single excitation states. Studying a system of two coupled cavities 
the authors of Ref. [13] introduced delocalized modes to reveal that the detuning parameters between the 
atoms and these modes govern the dynamics of the propagation. The possibility of control of the individual 
coupling constants between the adjacent cavities has been investigated in [14]. For the choice of a coupling 
that varies parabolically on the cavity sites a dispersion-free propagation of the single excitation states has 
been observed [2]. 

Following the procedure of [13 here we study the propagation of atomic and photonic states between 
two sites of an array of cavities with two different choices of coupling coefficients linking photons of adjacent 
cavities. Generalizing the results of |I3j we first consider a translation invariant closed chain of N identical 
cavities with periodic boundary condition where the coupling coefficient between the photons of nearest- 
neighbor cavities is held uniform. Introducing the delocalized atomic and photonic modes we follow a 
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time-averaging procedure [TSl [13] to perturbatively evaluate the effective Hamiltonian for the three hmiting 
cases: (i) dominant hopping parameter, (ii) large detuning between the atomic and photonic frequencies, 
and (iii) the case when the resonance between the atomic frequency and an eigenfrequency of a delocalized 
photonic mode is realized. It is interesting to note that in terms of delocalized atomic coordinates a diagonal 
spin-spin effective interaction that directly transfers energy without the intermediacy of a photonic process 
is introduced. 

Allowing for the precise control of the coupling constant between the adjacent sites the authors of Ref. 
[14] introduced the parabolic coupling scenario. In the context of one dimensional spin chains |T6j it has 
been observed [17] that the parabolic coupling between nearest-neighbor spins leads to exact transmission 
of a single excitation quantum state. In the limit of a dominant hopping parameter between the photons of 
the adjacent sites perfect transfer of a quantum state at pre-specified times is also realized in the present 
case. In the model studied here the atom-photon coupling at individual cavities introduces a dispersive effect 
that prohibits perfect transfer of a quantum state for large values of the said coupling. However, in another 
limiting case when the detuning parameter is large compared to both the hopping parameter and the atom- 
cavity photon coupling constant, exact propagation of the quantum state is also obtained with an increment 
in the required time. For these cases the exact time-dependent correlation function of the quantum state 
between two arbitrary cavities is obtained in a closed form. The plan of the paper is as follows. In Sec. |ll] 
we discuss the case of uniform coupling between the photons of adjacent cavities arranged periodically on 
a closed loop. Our discussion of the case of parabolic coupling between the photons of adjacent cavities is 
contained in Sec. IIIII Following this we conclude. 

II Cyclic cavities with uniform coupling 

Here we consider a cyclic chain of N identical cavities obeying periodic boundary condition, and each 
containing a two-level atom that is coupled to a localized photonic mode modeled as a harmonic oscillator. 
The transition frequency e of the atoms and the frequency of the oscillators are held uniform over the 
chain. The atom-cavity photon interaction is described by Jaynes-Cummings model with the adoption 
of rotating wave approximation. The adjacent cavities are interlinked via photon hopping with uniform 
couplings between neighboring cavities. The Hamiltonian may be expressed as a sum of the Jaynes-Cummings 
Hamiltonians for identical cavities and a nearest-neighbor photon hopping term: 

Af-l 

i/= F/c^^hop^ (2.1) 

j=o 

where the Jaynes-Cummings Hamiltonian for the j-th cavity and the hopping term between adjacent cavities, 
respectively, read: 

i?/*^ = e a+cr- + n a]aj + g (a+a^ + aja]^ , 

JV-l 

ff^°P = K ^ ajCj^fcafc, Cj^k ^ Sj k+i + Sj+ik- (2.2) 

In the above expression the iV-th and the 0-th degrees of freedom are identified. The standard commutation 
relations for the photonic and the atomic modes pertaining to individual cavities are given below: 

iaj,al] = ^jk, c^fc ] = ±2 6.jk erf, kj^, cr^ ] = S^k cr^k- (2-3) 

The operator Af that represents the total number of atomic and photonic excitations of the combined system 
commutes with the Hamiltonian (|2.ip : 

JV-l 

= E + ^h) ' [H,M]=0. (2.4) 

In this work we focus on the time-evolution of single excitation states that may be arbitrary superpositions 
of one excitation states of atoms and photons. The ground states of the photonic and atomic systems may 
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be listed as |0) = |{0j}), |G) = for j = 0, 1, . . . , iV— 1. The photonic and atomic one excitation states 

localized in j-th cavity are created as follows: 

a]|0) = |l,), a,|U.)-<5,,. |0), r7+|G) = |e,), aj\e,) ^ S, , \G) . (2.5) 

To diagonalize the photonic part of the Hamiltonian we need to introduce delocalized coordinates unitarily 
related to their local analogs pertaining to a cavity. For w being a root of unity the unitary transformation 
reads 

L/jfc = ^w^'' Vj,fc = 0, 1, a;^ = l oj ^ exp {2TTi/N) , (2.6) 



1 

- ^ t^J'C^-^) = 4f ^ UU^^U^U = I. (2.7) 



and the unitarity constraint readily follows 

N 



We note that here and elsewhere in this section we follow the modiV arithmetic. Now the delocalized 
photonic coordinates obeying canonical commutation relations are given by 

N-l N-l 

= 7^ E ^1 = 7^E^"''4 ^ [A,,Al]^5,,. (2.8) 

It turns out that the atomic coordinates that couple with the delocalized photonic coordinates in the recon- 
structed interaction Hamiltonian are also delocalized in nature. These collective spin variables are introduced 
as discrete Fourier transforms 

N-l ^ N-1 

Sf = ^ll^^''-t S^-^T.^^'^'^'k' iSfy=Sf, {S]y=Sl^ (2.9) 

^ k=0 k=0 

obeying a closed algebra: 

2 

N 

The unitary operator p.6p diagonalizes the coupling matrix between adjacent cavities given in (|2.2p : 



[S^, S^] = ± - S^^^ , [S+,S^] = S^_,. (2.10) 



(C/C C/t)^.^ = {lo^ + oj-^) 6,k = 2 cos (^j Sjk- (2.11) 
Recasting the Hamiltonian (|2.ip via the delocalized coordinates we obtain 

H^Ho + H^'''\ i7o - i/o^''^*^ + ^^0^^°*°"", (2.12) 

where 

N~l N-1 / o \ 

Hr''' = e S^S-, H^-'-- = n,A]A,, = n + 2. COS (, I ) . (2.13) 

j=0 j=0 ^ ^ 

For N > 2 the eigenfrequencies Vlj of the delocalized modes Aj are degenerate: ilj = VL^-j for j ^ 
{N — j) mod N. The non-degenerate eigenfrequencies correspond to the 'center-of-mass' mode Aq and 
'alternating' mode An/2 for even N. The delocalized coordinates introduced above maintains the structure 

that reflects the rotating wave approximation for the atom-cavity photon interaction term H^^^ in (|2.12p : 

N-1 

i/int=g ^ (^S+A,+S-Al). (2.14) 
The total number of excitation operator (|2.4p now assumes the form 

N-1 
3=0 
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The delocalized photonic and atomic single-excitation states may be introduced via invertible unitary 
transformation as follows: 

N-1 

^ 5+|G) = ^^c.-^'=|e,), S-\^k)=S,k\G), = 6,^- (2.16) 

^ k=Q 

The localized and the delocalized one excitation states introduced in p.5|) and (|2.16p are two sets of mutually 
unbiased bases: 

(l,|L) = -^c.-^■^ (e,.|^,) = ^c.-^'= ^ |(i^.|i,)| = |(e^.|^,)|= 1 . (2.17) 



The distinction between the localized and delocalized states may be summarized by introducing a unitary 
shift operator that translates localized photonic and atomic states through a single cavity. Specifically, it 
acts on the single excitation states as follows: 

r|l,) = r-i|l,)-|l,_i), r|e,) = r'^ |e,) = |e,_i) = r~\ (2.18) 

The construction (j2.16p makes it apparent that the delocalized states are invariant under the action of the 
shift operator r where the eigenvalues correspond to the root of unity phase angles: 

T|i,) r\^,)^io^i^,). (2.19) 

The delocalized photonic and atomic operators introduced in (|2.8p and p.9p transform as follows: 

tAj = uj-^ Aj, tA] = w^' A], t Sf t'^ = 5=^. (2.20) 

It is observed that the interaction Hamiltonian (j2.14p is invariant under the shift transformation: 

^^int^-i ^^int^ (2.21) 

The most general one excitation state may be equivalently expressed as linear compositions of either the 
localized or the delocalized basis states. These alternate expansions are listed, respectively, below: 

N-l 

\^{t)) = ^ (a,(<) |G) ® |1,) + b,(t) |e,) ® |0)) , (2.22) 

3=0 
N-1 

Invertible unitary transformations interrelate the coefficients of the above expansions (|2.22p and (|2.23p : 

N-1 N-1 

W = ^ E ^'^ W' W = ^ E b,(0, (2.24) 



N-1 , N-1 



The action of the Hamiltonian (|2.12p on the two dimensional subspace {|G) (E) <E) |0)} reads 

F|G)0|i,) - %|G)® |i,)+g|z9,)® |0), 

i?|i?,)®|0) = e|^,)®|0)+g|G)®|i,). (2.26) 
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The Schrodinger equation for the state {"^{t)) now block diagonahzes the one excitation subspace of the 
Hamihonian (|2.12l) in 2 x 2 blocks with the j-th block given by 



(2.27) 



The eigenvalues of the matrix Hj and the corresponding eigenvectors in the two dimensional subspace 
{|G) (g) <8) |0)} read, respectively, as 



%,± = — ^ — ±Xj, Xj = 



n, 



2k cos 



(j^), S = e-n, (2.28) 



3. ± - £J 



((%±-e) |G)®|i,)+g|^,)®|0) 



(2.29) 



The detuning parameter 5 measures the difference of the atomic frequency and that of the cavity pho- 
ton, whereas the corresponding differences for the delocalized modes are given by the parameters {Aj|j = 
0,1,...,7V — 1}. The time evolution of the coefficients in the expansion (I2.23P of the state |^'(i)) in the 
delocalized basis is given by 



aj{t) 



exp {—iHjt) 



a,(0) 



where the explicit construction reads 



exp 



■t 



P,{t) = cxp(-z^^t 



cos(xji) + « V 2 
cos(xj t) " M ^ 



. / A A sin(xji) 



X] 



sin (y it) 

Xj 

Xj 



(2.30) 



(2.31) 



To derive the effective Hamiltonian that discards the rapidly oscillating high-frequency components on 
account of a time-averaging process, we use the interaction-picture scenario. The interaction Hamiltonian 
given in (|2.12l) now assumes the form: 



(2.32) 



Its explicit evaluation following from (j2.12j|2.1"i|) reads 



JV-l 



H^'^^it) ^gJ2 ('=^P(* ^) ^t^^ + exp( - i t) 57 



(2.33) 



As we are interested in dynamical processes occurring at low frequencies, all high frequency components of 
the Hamiltonian are assumed to average out to zero resulting in a 'coarse-grained' effective Hamiltonian. 
Following the recipe given in [TSJ|T3] the time-averaged effective Hamiltonian up to the order 0{g^) reads: 



= Ho + (h"^*(<)) + i (\n^''\t),v{t) 



(2.34) 



where the averaged value of a dynamical quantity 0{t) with a suitable probability density function f{t) 
is defined [HI [l3] as (0(t)) = JZc>fi^ " t')0{t')dt'. The generating element V(t) of the order 0(g2) 
contributions in the rhs of (j2.34p may be obtained d la [15l [13] : 



N-l 



(2.35) 



3=0 



We now study various limiting cases of the effective Hamiltonian (|2.34p . 

The large hopping limit is dominated by the coupling between the photons of adjacent cavities: S,g. 
We also assume A^- » g V j € {0, 1, . . . , A^ — 1}. These limits ensure vanishing of the following time-averaged 
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quantities: (Ti^'^^{t)J = 0, (V(t)) = 0. Eliminating the high frequency components via the averaging 
process the commutator may be constructed using the definitions (I2.33P and (|2.35p : 



7{i"*(t),V(t) 



N>j>0 J 



SfAj,S 



(2.36) 



The second term in the rhs of (|2.36p owes its origin to the degeneracy of the eigenfrequencies that is present 
for > 2, and has been referred to in the context of (|2.13p . Exphcit evaluation of the commutators following 
from dm) and (|2J0)) 



[S+A, , S^_^A^^_^] = Sl^_^ A,A^j,_^ J^N-j mod N 
allows us to obtain the effective Hamiltonian up to the order O(g^): 



(2.37) 



eff 



E 

l->3>a 



N-1 

j=0 J 
1 

A" 



'0 {4 



A, 



(s+s,- + s,-s;)i 

)■ 



(2.38) 



The first term in the rhs of (|2.38p that may be read from (|2.12l) describes the energy of the bare delocalized 
atomic and field modes. The first sum represents the contribution due the Stark shift, and an induced 
interaction of order O(g^) operating directly between the delocalized spin modes without the intermediacy 
of the photonic modes. It is interesting to note that the delocalized modes maintain the diagonal structure of 
the spin-spin interaction term in the effective Hamiltonian. The second sum is present due to the degeneracy 
of the field modes. The expression (|2.38l) of the Hamiltonian via the delocalized modes makes it evident that 
there is no transfer of energy between the atomic and the photonic modes. Employing the reduction of the 
'averaged' spin coordinate 5*0 on the projectors Sq = [Pe — Pg), where P^ = J^J^o^ 



Xj >< Xj\yx e {g,e}, 



and with a suitable rearrangement of the terms in the first sum in (j2.38p we recast H'^^ as follows: 



eff 



Ho 



N 



P. 



N-l 



^ A, 



— A J a] Pc 



A\A, 



7V-1 
J=0 



■? k.i=0 



f >j>0 



N-l 



-2jfe A . A'\ 
-ti- n 



N-J 



-0.2^"'= A]Am-j 



k=0 



(2.39) 



In the above expression the Stark shift depends upon the population of the delocalized photonic modes. The 
coupling of spins between any two arbitrary distinct site exist at the order O(g^). For N — 2 the effective 
Hamiltonian (f239l) agrees with the result obtained in |13j at the corresponding limit. To make the process 
transparent we introduce a discrete Fourier transform on the reciprocal of the detuning parameters of the 
delocalized photonic modes: 



1 1 

G{j) = - y ^ 

k=0 



and thereby reexpress (|2.39|) in the following form: 



(2.40) 



eff 



N 



(N-l N-l 
Pe J2 ajGij-k)al~Pg ^a]G(j-fc) 
j,k=0 - '•-n 



ak + N 



j,k=0 



E 

f >j>0 



1 \ 

^3 k=o ^ ^ > 



N-l 

E". 



(2.41) 
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It is evident from the structure of (I2.4ip that the Fourier transform G(j) acts as the propagator of the 
atomic and the photonic excitations. The transition amphtudes of these excitations proceeding on the loop 
are proportional to the magnitude of the propagator. 

To further examine the issue of transfer of the excitation modes we turn to the solutions (|2.3ip . The 
following limits valid in the present approximation 



A, 
2X, 



g<- 



Xi 



g<- 



(2.42) 



yields the structure 

a,{t) = exp (^-z (n, - ?j a,iO), (3,{t) - cxp (^-z (^s + ^) ?j /3,(0). (2.43) 

It follows from the above equations for the delocalized modes that the atomic and photonic excitations 
remain decoupled from each other. Employing these solutions with the Fourier expansions (|2.24[ I2.25P we, 
in this limit, obtain the time-evolution equation for the coefficients of the localized single excitation states: 



I. c—n \ \ 



kl=0 
N-l 



k.f.=t) ^ ^ 



m 



afc(O), 
bfe(O). 



(2.44) 



Equation (j2.44p indicates that subject to the approximation (I2.42p transfer of excitations from the atomic 
to the photonic modes or vice versa does not take place, and, consequently, the probability density in each 
kind of degree of freedom is conserved: 



N-l N-l N-l N-l 

|a,(0|^ = E |a,(0)^ Ib.Wl' = E Ib.(0)l' 

]=0 j=0 j=0 i=0 



(2.45) 



For definiteness, we assume that at t = only the 0-th atom is in the excited state and all other degrees of 
freedom are in the ground state: aj(0) = 0, bj(0) = (5jo. Evolution equations in (|2.44p now reduce to 



1 r 

a,(0=0, b, W = -^ cxp 

fc=0 ^ 



— i I jk — 



(2.46) 



For N ^2 exact transmission of the excitation does not take place. 

We next study the large detuning limit: (5 3> k, g when the atoms are highly detuned from the photonic 

modes. In this limit also the following time-averag ed properties hold: (n^^^{t)^ = 0, {V{t)) = 0. Towards 
obtaining the effective Hamiltonian we first evaluate the time-averaged commutator in this limit: 



H^n*(t),V(t) 



N-l 



N-l 



3=0 



j,fc = 



1 



1 



Afc 



(2.47) 



Employing the first equation in (|2.37p and rearranging terms on the rhs of (|2.47p the said commutator is 
reexpressed as 



( N-l 



N-l 



1 



■ ■ " ' i=o ^ 



j",fe=0 ^ 

The effective Hamiltonian up to the order O(g^) may now be obtained via (|2.34l) and (I2.48P : 



(2.48) 



( N-l N-l 

E TT [^U A A + A]A^ + Yir (St Si + SJS^) I 



— A, 



A,- 



(2.49) 
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In the N = 2 case the above expression reduces to the results obtained in [TJ] in the large detuning limit. 
The effective Hamiltonian (|2.49l) may be recast using the propagator in the Fourier space (|2.40p . and the 
transforms 



n 



eff 



2 JV — 1 y 

^" + Y E {a]GU-k)ak + a,GU~k)4) 

+ (a+ G{j - k) aZ + a-G{j~k)a+) l] . 



(2.50) 



In this approximation also the atomic and the photonic modes remain decoupled from each other forbidding 
any transfer of energy up to the order O(g^). The time-evolution of the coefficients aj{t),hj{t) of the 
localized single-excitation states reduces to previously described form (|2.44l) . Assuming dominant value of 
the detuning parameter 5 ^ k, g we now study the asymptotic N ^ oo limit of p.44p for the boundary 
condition aj(0) = Sjo, bj{0) = 0. Using the Fourier transform 



exp(±ia; cos 6*) — (i*)" Jnix) exp(±w 



Jn{x) = i-ir Jn{x) 



(2.51) 



via the Bessel functions Jn{x), and employing the said boundary condition the time evolution of the coeffi- 
cients of the localized states may be expressed as 



(2 \ ^ { 2 \ 



W = 0. 



(2.52) 



In deriving (|2.52p we have used the summation (|2.7p . The asymptotic expansion of the Bessel function of 
large order [18] Jn[]^)\n^oo ^ N~^/^ allows us to obtain an asymptotic scaling limit of the magnitude of 
the transition coefficient |a^/2(^ ^ Af)|Ar-i.oo ~ N~'^/'^. 

The case of comparable hopping and detuning is of particular interest. Turning towards this near- 
coherence scenario where the atomic frequency approximately equals to the eigenfrequncy of, say the ^-th 
delocalized photonic mode: e ~ =^ 0, we observe that the degeneracy (|2.13p of the photonic 

eigenmodes present when I ^ N — I mod needs to be taken into account. We also assume that the 
detuning parameters for all the nonresonant modes are large: ^ Vj ^ {£, N — £ mod N)}. Unlike 
the previous cases considered the present near-resonance limit allows a direct transfer of the excitation 
between the atomic and the photonic modes. The time-averaged contribution of the near-resonance mode 
to the T-O^^ at the order 0(g) may be read from ()2.33p : 



H"^*(t) 



,N-i thoAn 



) {"^N-iAj^^i + S 



4' 



)) 



whereas its contribution to (V(t)) follows from (|2.35p 

{m)^^[{siAi 



St A, 



(l '^^,Ar-£ modw) {^N~i^^N-l 



S 



N-iAN-i 



(2.53) 



(2.54) 



When a near-resonance scenario involving these degenerate eigenmodes are realized the said modes give rise 
to independent contributions associated with them, respectively. The commutator of the above time-averaged 
quantities reads 



y}''\t)),{V{t)) 



2g^ 



St Ai^ So A\ 



+ (1 - 5 



mod Af) ^ 



A 



NJ 



S^i-eAN^i, Sj^_gA]^_g 



~^ S2£-]^ A(>A^j^_g Sj^_2i A^gA]\i—g. (2.55) 



Combining the previous results the effective Hamiltonian (|2.34p in the near-resonance case up to the order 
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O(g^) assumes the form 



eff 



Ho+g[ {S+A, + Si 4) + (1 - 6,^^_, ^^odiv) (S^-i^N-^ 

2 

+y E' ^ (^0 {aIa, + A, A]) + {s;s- + s-s;) l) 

j=Q ^ 



+ 'S'jY_^Ajy_^) j 



(2.56) 



->j>o 



where the first sum on the rhs excludes the resonance modes {£, {N — £) when £ ^ {N — £) mod N}, and the 
resonance mode {£) lying in the domain of the index of the second sum is eliminated. The term 0(g) in this 
Hamiltonian describes a transfer of energy between the near- resonant delocalized photon modes and their spin 
excitation partners. The non- resonant modes contributing in the first sum at order O(g^) represent the Stark 
shift of the atoms and a direct transfer of excitations between the atomic modes without the intermediacy of 
the field modes. In the near-resonance condition we assume the detuning parameter of the resonant modes 
are negligible compared to the atom-cavity photon coupling: ^ <C 1 Vj G {£, {N — £) when £ (N — £) 
mod N}. Excitations of the near- resonant field modes and that of their coupled spin partners may be read 
from p3T|) : 



,{t) = cxp -qe- 



A, 



I3j{t) = exp(-i(e-^)t 



cos H g + ^) *) " * ^^'^ ((g + ^ ) M ^^(0) 

cos 1 1 g + ^) - * ((g + ^ ) M "^(0) 



(2.57) 



Assuming that the detuning parameters for the non-resonant modes satisfy 



A, 



£) when £ ^ {N — £) mod N}, the corresponding excitations assume the form 



< 1, where j ^ {£, {N 



The above evolution equations for the delocalized modes in conjunction with the Fourier transforms (j2.241 
I2.25P now yield the single-excitation states of the atoms and the cavity photons. We first assume that the 
near resonant £-th. mode is nondegenerate: £ = {N — £) mod A^. The i = 0, N/2 (for even TV) modes satisfy 
this property. For the choice of the initial excitations being atomic in nature: 3^(0) = 0, the time evolution 
of the coefficients are given by 



b,(i) = 



exp 



N 

\11 E 



-I e 



2 



sin I 



N 



k=0 



a = 



-l-w"^^-''"''^ exp 



'I e 



At 
2 



'I e 



t cos 



8i 



^ k=0 



A2 

8g 



)')] 



bfc(O). 



(2.58) 



In the exact resonance condition = 0, and in a scenario where the dispersive effects are negligible 



A,- 



y j ^ £ the above solutions assume the form 



N-l 



3 At) = -J7 cxpi-zet) sin(gt) ^ u;-'^^-'^^ bfe(O), 



A^ 

exp(— zei) 



k=0 



b,(0)-^sin^(|) E--'^^-'^b..(0) 



(2.59) 



With the further assumption that only the 0-th atom is initially excited whereas all other atoms remain in 
the ground state bj(0) = Sjo Vj G {0, 1, . . . , A^ — 1}, the coefficients of the localized single-excitation states 
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read 



ajW = -— sin(gi) 



N 

exp(— iet) 6jo 



cxp I -i (j^'^ + j 



(2.60) 



For N ^ 2 the magnitude of the locahzed atomic excitation bi(<) assumes unit value for both choices of the 
resonance modes ^ = 0, 1 at the specified time t = {2n + l)7r/g Vn G Z+, making the transmission of the 
state exact. 

Turning towards the resonance of the atomic mode with a conjugate pair of degenerate delocahzed 
photonic modes {£, {N — £)}, where £ ^ (N — £) mod iV, we, following the preceding recipe, construct the 
time-evolutions of the excitations of the atoms and the cavity-photons: 



2i 



-i \ s - ^ \ t\ sin 



fc=0 



^-"(j-fc) gxp 



= 



+2 exp 



-I e 



2 



t cos 



bfc(O). 



(2.61) 



In the above equation we have assumed that the initial excitations are only atomic in nature. For the exact 
resonance condition Ai = and in the dispersion- free limit ~ \f j ^ {£, {N ~ £)} the excitations reduce 
to the form 



aj{t) = cxp{-iet) sin(gf) ( 7^ ~ ) ^f"^^^ 



N-l 



b,(t) 



N 

exp{—iet) 



k=0 



2tt 



N 



b,(o) - - sin^ ( ^ ) E con ^ - ^) ) ^'^^^ 



k=0 



2n 



N 



(2.62) 



Restraining to the case where only the O-th atom is initially excited bj(0) = Sjq Vj G {0, 1, . . . , — 1} the 
time evolution of the coefficients read 



sin(g<) I exp (^i (^j£^ - et 



cxp ( -i ( ^ + et 



bj W = exp(-iet) (5jo - 



cxp ( i ( - et 



cxp ( ( + £t 



(2.63) 



where a superposition of the clockwise and the anticlockwise modes takes place. For = 4 and with 
the choice of the resonance mode £ — 1 the magnitude of the excitation h2{t) assumes unit value at i = 
(2n-h l)7r/g Vn G Z+. 



Ill Linear chain with paraboUc coupUng 

Interesting situations arise if one considers nearest-neighbor couplings between the cavities to be nonuniform 
in nature. In particular it has been observed P3„ that a parabolic coupling between the photons of adjacent 
cavities gives rise to dispersion free transmission of the excitations along one dimensional chain of cavities. 
The inter-cavity photonic hopping term of the Hamiltonian now reads 

N-l 

i/hoP = « 51 «] ^J.fe "-t' = VjiN-j) S, k+i + VHN-k) S,+i k. (3.1) 

The photonic part of the interaction Hamiltonian relates to the tridiagonal Jacobi matrix of the Krawtchouk 
polynomials [19]. Therefore the photonic degrees of freedom may be diagonalized using delocalized wave 
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functions expressed via these discrete orthogonal polynomials. Very briefly we now introduce the standard 
notations on the Krawtchouk polynomials. Detailed discussions may be obtained from Ref. 

The Krawtchouk polynomial of degree n {n — 0,1, . . . ,W) in the variable x, with parameter < p < 1 is 
given by 

Kn{x) = Knix;p,m) = 2Fii (3.2) 

The function 2F1 is the classical hypergeometric series [20l [21], and in this case it is a terminating series 
because of the appearance of the negative integer (— n) as a numerator parameter. It is convenient to 
introduce orthonormal Krawtchouk polynomials by 



K„{x) ^ ^l^^EM, (3.3) 
where w{x) is the weight function in x, and d„ is a function depending on n: 

w{x)={^^p^{l-p)'^-^ Vx = 0,l,...,ai; d„ = -py(^i^) . (3.4) 

The scaled polynomials Kn{x) satisfy a discrete orthogonality relation [19 : 

^ Kn{x)Kjn{x) = Snm- (3.5) 

In the present scenario of parabolic coupling of the photons of adjacent cavities the Krawtchouk polyno- 
mials of parametric value p — 1/2 play an essential role in diagonalizing the photonic part of the Hamiltonian. 
Following [55] we introduce a symmetric orthogonal N x N matrix that is comprised of Krawtchouk poly- 
nomials as 



Ujk^Kk(^j;^;N-l^ V j, fc = 0, 1, . . . , TV - 1, U - U^, UU^ = U^U=I. 



(3.6) 



The adjacency matrix C given in p.ip may now be diagonalized via the orthogonal matrix U as follows: 

UCU^^D, D = diag{N -1, N -3,...,-{N -I)). (3.7) 

We note that the parity symmetry of the scaled Krawtchouk polynomials described below 

interrelates the components of the orthogonal matrix U: 

Ujv-i, = (-1)^ Uo,. (3.9) 

The delocalized collective modes may now be introduced as transform of the corresponding localized 
degrees of freedom via the orthogonal symmetric matrix U as given below: 



k=0 k=0 

N~l N-1 N-1 

S+ = Eu.fc^^ S7 = ^U,.c7,7, S^^,= Y.U,iUMa! ^ [S+, S,t] = S|,. (3.10) 



fc=o fe=o e=o 



The photonic degrees of freedom {Aj\j £ 0,1, ... N ~ 1} satisfy the Heisenberg algebra. We now recast 
the Hamiltonian employing these delocalized variables. The rotating wave structure of its atom-photon 
interaction term is maintained in terms of the delocalized modes: 



^0 = s E S; S7 + E %A] A„ H-t ^ g J2 (S; A, + ST At) , (3.11) 

3=0 j=0 j=0 
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where the spectrum of the delocaUzed photonic eigenmodes is given hy ilj ^ fl + k {N — 1 — 2j) for 
j = 0,1, . . . , N — 1. The delocahzed single excitation atomic and photonic states are constructed by the 
action of the collective operators introduced in p.lOp : 

JV-l 
fe=0 

= S+|G>= ^U,,|efc), Sj\h)=S,k\G), {d,m=S,,, (e,|4) = U,,. (3.12) 

k=0 

In the above delocahzed basis an arbitrary single excitation state may be expanded as follows: 

Af-l 

l*W> = E ("J-W ® + W 1^.) ® |o)) . (3.13) 

Invertible transformation via the orthogonal symmetric matrix U interrelate the coefHcients of the above 
expansions with the the corresponding coefficients of the expansion (j2.22p in the basis of the localized single 
excitation states: 



ij(<) = E Ujfeafe(t), /3,(t) = E Ujfebfc(i), a,(<) = ^ Ujfeafe(i), h,{t) = ^ ^ikPkit). (3.14) 



N-l JV-1 W-1 N-1 

k=0 fe=0 k=0 k=0 

We first consider the regime of dominant hopping parameter: S, whereas the atom-photon coupling 

is assumed to be small: ^1 Vj € (0, 1, . . . , — 1). Unlike the case considered in Sec. HIl the 

parabolic coupling between the photons of adjacent cavities does not produce degenerate eigenvalues of the 
diagonalized photonic modes. The effective Hamiltonian that eliminates the rapidly varying frequencies may 
be obtained following the recipe described in Sec. [TTl We quote the result below: 

^'^ = Ho + \Y.x- (^^^ (^'^^ + ^^"^O + ^^^^ + ^^^^^ ^) ' (^-'^^ 

where the detuning parameters for the photon eigenmodes read Aj = e — Vlj = S — k {N — 1 — 2j). Towards 
obtaining the transmission of single excitation quantum states we study the evolution of the coefHcients 
aj{t), bj{t) of the localized excitations. In the present limit the equations p.l4p in conjunction with (|2.3ip 
produce the time-dependent coefficients as 

Af-l N-l 

a,{t) = E afe(O)' bj(t) - E bfc(O), (3.16) 

k=0 k=0 

where the kernels encoding the correlation functions read 



JV-l X / 2 \ \ 

^a(j, k;t) = E ^ji Ufe£ exp - -^j tj 



ICbij, k;t) = E Ufe£ exp ie + |- j t j . (3.17) 



1=0 



Of particular significance is the case where the dispersive coupling between the atoms and the corresponding 
localized cavity excitations is neglected: g = 0. The evolution kernel of the localized photonic state reads as 



N-l 



/C3(j,fc;i)=exp (-i(r! + K(iV-l))t) /j-fc(i) /,,fe(t) = ^ 2 = exp(z2Kt). (3.18) 



1=0 



By employing a classical summation [23| formula of the hypergeometric series the time-dependent correlation 
between the j-th and k-ih. sites may be expressed in a closed form : 



N-\\fN-l\,^ -j,-k -Az 
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In the context of a linear spin chain with the interaction determined by a Jacobi matrix such propagator of 
a single excitation state was earlier obtained in [53]. Assuming that the initial state in p.l9p is given by 
a photon in 0-th cavity: 3^(0) = Sjo, the transmission of the photonic excitation to an arbitrary cavity is 
described by the evolution of the said coefficient as follows: 



a,(t) = i-^y 




exp(— ifit) (sin(Kt))"' (cos(Kt)) 



N-l-j 



(3.20) 



It is evident from the above structure that in the dispersion-free g = case the amplitude of the photonic 
excitation at the (N — l)-th cavity is of unit magnitude |aj(i)| = 6j n-i at the following predetermined 
times: {i = t„ = '■^"^^-'^ |n = 0, 1, . . .}. 

In the presence of the coupling between the atom and the cavity photon g 7^ it is not possible in general 
to obtain a closed form expression for the time-evolution of the coefficient of the single photonic excitation 
state. To demonstrate that an exact transmission for this case is not possible in the present limit k ^ 5, g 
we proceed as follows. For the choice of initial state aj(0) — 6jo the evolution kernels given in p.l6[ I3.17P 
may be employed to express the magnitude of the excitation in the (N — l)-th site as given below: 



\^N-l{t)\ 



N-1 



^ UN-i^Hoi exp [i(2«:^+l-)i) 



(3.21) 



The parity relation (13. 9p allows us to recast p.2ip as 



|aw-i(OI 



N-l 2 ( 2 \ 

^ (Uof) exp z[(2n+ l)7r£-H(2K£ +!-)<] 
i=o ^ J 



0,1,.. 



(3.22) 



The relative phase factors of the terms in the above summand do not, in general, identically reduce to zero 
for any choice of time. Consequently, the magnitude of the excitation aN-i{t) obeys the inequality 



|ajv-i(i)l < 



Af-l 



E (uo.)' 



= 1. 



(3.23) 



Only in the vanishing limit of the coupling g between the atom and the cavity photon, and also for 
large detuning limit (5 3> K,g to be discussed next the excitation |aAr_i(<)| equals its unit limiting value for 
specific predetermined values of time. 

Proceeding parallel to the derivation in Sec. |ll] the effective Hamiltonian in the limit of large detuning 
parameter S ^ K,g may be obtained. Up to the order O(g^) the effective Hamiltonian reads 



n 



eff . 



JV-l 



j,k=0 



n A, 

3=0 1 



's+s- 

y J 3 



3 3 



(3.24) 



The propagation of single excitation states follows the description given in (|3.161 13.17p . To evaluate the 
propagator we expand the exponent in (I3.17P retaining terms up to the order 0{j), and neglecting higher 
order terms in the said coefficient. The kernel A^a(j, /c;t) given in (j3.17p may now be evaluated exactly as 
before reproducing the result p. 18) I3.19P with a redefinition of the constants: 



K 1 



We reproduce the kernel /Ca(j, fc; t) as follows: 

/Ca(j,fc;i) =exp ^-i(r2' + K'(7V-l))ij ^\JjiUki{z'Y, z' ^ cxj>{i2K't). 



N-l 



(3.25) 



(3.26) 



Therefore, as in the dispersion free limit discussed in p.20p . the exact transmission of a single photon 
state from the 0-th to {N — l)-th cavity also takes place when in the exponent of the relative phases in 
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the correlation function p.l7p we retain only the linear terms in the parameter j. The time interval of 

exact propagation of the state now assumes the value = ^^""^l^^ Vn = 0, 1, . . .. In the context of the 
approximation used here the effect of a nonvanishing value of the coupling constant g 7^ is to increase the 
time of transmission of the single excitation state i. e. to decrease the velocity of propagation of the quantum 
state. Retaining terms of quadratic or higher order 0((f)^) in evaluating the kernel IC3{j,k;t) defined in 
p.l7p we observe that the exponent depends on the summation index at least quadratically. Consequently, 
the summation no longer remains of hypergeometric type, and therefore a general closed form expression is 
not obtained by the present technique. 

Lastly we study the resonance limit e « fi^ « for the £-th delocalized photon eigenstate in the 

presence of the parabolic coupling between the photons of adjacent cavities. We also simultaneously assume 



that the following hierarchy of values holds: 



^ 1 Vj 7^ In contrast to our description of 



« 1, 

the resonance scenario given in Sec. [Tl] the delocalized photon eigenstates here are not degenerate. In the 
resonance limit the effective Hamiltonian has contributions of order 0(g) that couples atomic excitations 
with the photonic excitations: 

H^ff = Ho+g f S+ A, + S7 A| 



N~l 



Ej-i^h (A.A] + A]A,) + (S;S7 + S-S+)l) . (3.27) 



2 U A, 

Assuming that the initial excitation is only atomic in nature: 3^(0) = Vj S (0, 1, . . . , iV — 1) the time- 
evolution of the coefficients may be obtained via p. 141 12.3ip as given below 

a,(i) = -i exp (^e- ij sin|(^g+^^tj ^ Uj£Ufc^bfc(0), 



bj(*) = X! X! Uj„Ufc„exp 

^— n ^ Ti=o 



fc=0 "=o 



+ Uj, Ufc, exp ( ( £ - '"'^^^^^ ^^'^^^ 



For the exact resonance case = 0, and in the absence of dispersive effects ^ ^ V j 7^ £ the above 
evolution equations with the initial condition bj(0) = <5jo read 

aj (t) = -iUo£Uj^exp(-iei) sin(gi), bj(i) = exp(-iet) j,5jo - 2 Uq^ U^^ sin^ (^^^ j . (3.29) 

In particular the excitation of the {N — l)-th atom reads 



hN-i{t) = -2 exp(-zei) Uq, \^n-ii sin^ ( ^ ) = 2 (-1)^+^ exp(-iei) (Uq^)" sin^ ( |- ) , (3.30) 



where in the second equality we have used the parity relation p.9p . As it may be observed from the 
construction p. 61 13. 3p of the orthogonal matrix U that in the exact resonance regime for the parabolic 
coupling a perfect transmission of the quantum state is realized for the length of the array iV = 3, and for 
the resonance mode ^ = 1 at the time (^"+^)'^ Vn G (0, 1, . . .). 



IV Conclusion 

Here we have considered one dimensional array of optical cavities with a nearest-neighbor hopping interaction 
of the cavity photons described by Jaynes-Cummings-Hubbard model. In particular we have studied the time 
evolution of single excitation states in such chains of coupled cavities as models for quantum communication. 
Two different choices of coupling coefficients linking photons of adjacent cavities have been investigated. 
Employing delocalized collective photonic and atomic modes we first study a translation invariant closed 
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chain of an arbitrary number of N identical cavities with uniform couphng constant. Fohowing this we have 
considered a hnear chain of cavities where the non-uniform parabohc hopping term of the Hamiltonian is 
related to the tridiagonal Jacobi matrix associated with the Krawtchouk polynomials. For both of these cases 
we obtain the effective Hamiltonian in various simplifying limits. For dominant value of photonic hopping 
parameter, as well as for large detuning parameter between the atomic and the photonic frequencies the 
excitations are transferred between, say, the atoms without populating the field modes. Expressed via the 
delocalized atomic modes the spin-spin interaction term in the effective Hamiltonian is diagonalized. The 
atoms experience a Stark shift dependent upon the population of the field modes. In the case of resonance 
between the atoms and a delocalized photonic mode propagation of excitation between the atoms requires 
intermediate excitation of a photonic mode. In the case of non-uniform parabolic coupling between the 
photons of adjacent cavities the transmission of single-excitation states is exact for the limiting value g — 
of the atom-cavity photon coupling constant. However in the large detuning limit where we retain only 
the linear terms in the parameter j the propagation of the single excitation states remains exact with an 
increase in the time of transmission. This is observed in our evaluation of the time-dependent correlation 
function of the one-excitation states. 

The analysis developed here may have applications in certain contexts. The idea of having optical 
lattice systems with pre-engineered coupling constants between individual lattice sites may give rise to 
interesting physical situations. In particular the resonance situation where mixing of atomic and field modes 
is realized merits attention. In the present model it may be possible to generate multipartite entangled 
states of polaritonic qubits and study the time variation in the entanglement in the presence of dissipative 
atom-cavity photon coupling. Moreover, in the context of one dimensional spin chains it has recently been 
observed [25' that coupling between lattice sites determined by the Jacobi matrices of q-deformed Krawtchouk 
polynomials ensure perfect transmission of quantum states. It should be useful to study the problem in the 
context of Jaynes-Cummings-Hubbard lattices where coupled atomic and photonic modes are considered. 
The delocalized collective modes discussed here may allow easy extraction of the effective Hamiltonian at the 
desired perturbative level. Another possible application of the present method may lie in the quantum phase 
transitions in coupled array of atom-cavity photon systems |26j . Specific pre-engineered coupling between 
the cavity sites may give rise to desirable properties for a quantum transition between the Mott insulator 
and the superfluid states. 
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